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Abstract 

We have studied the color dipole picture for the description of the deep inelastic 
process, mainly the structure functions which are driven directly by the gluon distri- 
bution. Estimates for those functions are obtained using the Glauber- Mueller dipole 
cross section in QCD, encoding the corrections due to the unitarity effects which 
are associated with the saturation phenomenon. Frame invariance is verified in the 
calculations of the observables when analysing the experimental data. 
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1 Introduction 



In the kinematical region of small proton momentum fraction x, the gluon 
is the main parton driving the behavior of the deep inelastic quantities. The 
standard QCD evolution [1] furnishes a powerlike growth for the gluon distri- 
bution. This result leads, at first glance, to the unitarity violation at asymp- 
totic energies, requiring a sort of control. In the partonic language, at the 
infinite momentum frame, the small x region corresponds to the high par- 
ton density. The latter is connected with the black disk limit of the proton 
target and with the parton recombination phenomenon. These issues can be 
addressed through a non-linear dynamics beyond the usual DGLAP formal- 
ism. The complete knowledge about the non-linear dynamical regime plays an 
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important role in the theoretical description of the reactions in the forthcom- 
ing experiments RHIC and LHC, where these effects will be enhanced by the 
high energies or by the nuclear probes. 

The description of DIS in the color dipole picture is somewhat intuitive, pro- 
viding with a simple representation in contrast to the involved one from the 
Breit (infinite momentum) frame. Considering small values of the Bjorken 
variable x, the virtual photon fluctuates into a qq pair (dipole) with fixed 
transverse separation r at large distances upstream of the target, and inter- 
acts in a short time with the proton. More complicated configurations should 
be considered for larger transverse size systems, for instance the photon Fock 
state qq + gluon. An immediate consequence from the lifetime of the pair 
(l c = l/2m p x) to be bigger than the interaction one is the factorization be- 
tween the photon wavefunction and the cross section dipole-proton in the 7* p 
total cross section. The wavefunctions are perturbatively calculable, namely 
through QED for the qq configuration [2] and from QCD for qqG [3]. The 
effective dipole cross section should be modeled and it includes perturbative 
and non-perturbative content. However, since the interaction strength relies 
only on the configuration of the interacting system the dipole cross section 
turns out to be universal and may be employed in a wide variety of small x 
processes. 

We have taken into account a sound formalism providing the unitarity correc- 
tions to the DIS at small x, namely the Glauber-Mueller approach in QCD. 
It was introduced by A. Mueller [4], who developed the Glauber formalism to 
study saturation effects in the quark and gluon distributions in the nucleus 
considering the heavy onium scattering. Later developments obtained an evo- 
lution equation taking into account the unitarity corrections (perturbative 
shadowing), generating a non-linear dynamics which is connected with higher 
twist contributions. Its main characteristic is to provide a theoretical frame- 
work for the saturation effects, relying on the multiscattering of the pQCD 
Pomeron. In this contribution we report our studies considering the parton sat- 
uration formalism to describe the observables driven by the gluonic content of 
the proton in the color dipole picture [5]. The inclusive structure function F 2 
is calculated, disregarding the approximations commonly considered in pre- 
vious calculations [6]. The structure functions F L and F| c are also presented 
using the Glauber-Mueller approach and the rest frame in comparison with 
the experimental data. 



2 The DIS in the Rest Frame and the Glauber-Mueller Approach 

The rest frame physical picture is advantageous since the lifetime of the photon 
fluctuation and the interaction process are well defined [7]. The more simple 
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case is the quark- ant iquark state (color dipole), which is the leading config- 
uration for small transverse size systems. The well known coherence length 
is expressed as l c = l/(2xm p ), where x is the Bjorken variable and m p the 
proton mass. An important consequence of this formulation is that the pho- 
toabsortion cross section can be derived from the expectation value of the 
interaction cross section for the multiparticle Fock states of the virtual pho- 
ton weighted by the light-cone wave functions of these states [2]. That cross 
section can be cast in the quantum mechanical factorized form, 

1 

a£l(x, Q 2 ) =Jd 2 rJ dz \^ L (z, r)| 2 a di ^(x, z, r) , (1) 

o 

The formulation above is valid even beyond perturbation theory, since it is 
determined from the space-time structure of the process. The ^t,l(z, r) are 
the photon wavefunctions (for transverse T, and longitudinal L, polarizations) 
describing the pair configuration: z and 1 — z are the fractions of the photon's 
light-cone momentum carried by the quark and antiquark of the pair, respec- 
tively. The transverse separation of the pair is r. The explicit expressions for 
the wavefunctions are well known, 

\* T (,,r)| 2 = ^ £ e\ { [z 2 + (1 - zf]e 2 Kl (er) + m\ K 2 (er) } 2 ) 
\*l (*, r )! 2 = E e 2 {AQ 2 z 2 (l- zf K 2 (er)}. (3) 

i 

The auxiliary variable e 2 = z(l — z)Q 2 +m 2 , with m q the light quark mass, and 
Kq and K\ are the Mc Donald functions of rank zero and one, respectively. 
The quantity o- dl P° le is interpreted as the cross section of the scattering of the 
effective dipole with fixed tranverse separation r [2]. The most important fea- 
ture of the dipole cross section is its universal character, namely it depends 
only on the transverse separation r of the color dipole. The dependence on the 
external probe particle, i.e., the photon virtuality, relies in the wavefunctions. 
In general, an ansatz for the effective dipole cross section is obtained and the 
process is analized in the impact parameter space. The main feature of the 
current models in the literature is to interpolate the physical regions of small 
transverse separations (perturbative QCD picture) and the large ones (Regge- 
soft picture). Here we have used the Glauber-Mueller approach to determine 
the dipole cross section, with the advantage of providing the corrections re- 
quired by unitarity in an eikonal expansion. For the large r region, we choose 
to follow a similar procedure from the saturation model (GBW) [8], namely 
saturating the dipole cross section (r-independent constant value). 
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Now, we shortly present the main results from the Glauber-Mueller approach. 
Considering the scattering amplitude dependent on the usual Mandelstan vari- 
ables s and t, now written in the impact parameter representation b, 



«s,b) = ± J d 2 qe-^- b A(s,t = -q 2 ). (4) 

the corresponding total and elastic cross sections (from Optical theorem) are 
rewritten in the impact parameter representation (b) as 



atot = 2 / d 2 h Im a(s, b) ; a d = J d 2 b \a(s, b)| 2 , (5) 

An important property when treating the scattering in the impact parameter 
space is the simple definition for the unitarity constraint [6]. If the real part 
of the scattering amplitude vanishes at the high energy limit, corresponding 
to small x values, the solution to the that constraint is 



a(s, b) 



1 -e' 



fi(s,b) 



O'tot 



d 2 b 



l-e 



fi(s,b) 



(6) 



where the opacity Q is an arbitrary real function and it should be determined 
by a detailed model for the interaction. The opacity function has a simple 
physical interpretation, namely e~ n corresponds to the probability that no 
inelastic scatterings with the target occur. To provide the connection with the 
Glauber formalism, the opacity function can be written in the factorized form 
f2(s,b) = il(s)S(h), considering S(b) normalized as / d 2 hS(h) = 1 (for a 
detailed discussion, see i.e. [9]). 



We identify the opacity Q(s ~ Q 2 /x; r) = er 



nuclcon i 



x, r). The (qq pair) dipole- 



proton cross section is well known [6,9] and in double logarithmic approxima- 
tion (DLA) has the following form 



<^ m ^r) = ^f^r 2 xG{x^Q 2 ) (7) 



with the r-dependent scale Q 2 = r^/r 2 . Considering Eq. (7) one can connect 
directly the dipole picture with the usual parton distributions (gluon), since 
they are solutions of the DGLAP equations. In our case, we follow the cal- 
culations of Ref. [6,9] and consider the effective scale Q 2 = 4/r 2 . From the 
above expression, we obtain a dipole cross section satisfying the unitarity con- 
straint and a framework to study the unitarity effects (saturation) in the gluon 
DGLAP distribution function. Hence, hereafter we use the Glauber-Mueller 
dipole cross section given by 
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<ole = 2 Jd% (l-e^^o n (x,r)S( b )) 



(8) 



In order to perform numerical estimates one needs to define the profile func- 
tion S(b). This function contains information about the angular distribution in 
the scattering. We have chosen a simple gaussian shape in the impact param- 
eter space, S(b) = -^re~ b / r a, where A is the atomic number and Ra is the 
target radius. We will keep this notation although we are only concerned with 
the nucleon case. The R\ value should be determined from the data, ranging 
between 5 — 10 GeV~ 2 for the proton case [6]. Here, we have used the value ( 
R\ = 5 GeV~ 2 ) obtained from a good description of both inclusive structure 
function and its derivative [10]. Such a value corresponds to significative uni- 
tarity corrections to the standard DGLAP input even in the current HERA 
kinematics. 

In the calculations we have used the GRV94 parametrization [11]: bearing in 
mind that Q 2 = 4/r 2 , its evolution initial scale = 0.4 GeV 2 allows to scan 
dipole sizes up to r cut = ^ GeV -1 (= 0.62 fm). For recent parametrizations, 
where Qq ~ 1 GeV 2 (r cut ps 0.4 fm), the uncertainty due to nonperturbative 
content in the calculations would increase. An additional advantage is that 
GRV94 does not include non-linear effects to the DGLAP evolution since it 
was obtained from rather large x values, i.e. this ensures that GRV94 does not 
include unitarity corrections in the initial scale. To proceed, for the large r 
region, we choose the following ansatz: the gluon distribution is frozen at scale 
r cut , namely xG(x, Q 2 ut ). Then, for the large distance contribution r < r cut 
the gluon distribution reads as xG(x,Q 2 < Ql) = Q 2 /QlxG(x,Q 2 = Ql) 
leading to the correct behavior xG(x, Q 2 ) ~ Q 2 as Q 2 — > and providing the 
simplest technical procedure. 



3 Obtaining the Structure Functions 

3. 1 The structure function F 2 

First, we perform estimates for the structure function F 2 at the rest frame 
considering the Glauber- Mueller dipole cross section [5] . The expression, with 
the explicit integration limits on photon momentum fraction z and transverse 
separation r is, 

(l 2 °° 1 

H^Q 2 ) =T^— f d2r l dz (l*T(*,r)| 2 + |* L (z,r)| 2 ) ^f olc (x,r 2 ). 
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Fig. 1. The Glauber-Mueller (GM) result for the F2(x,Q 2 ) structure function. It is 
shown the transverse contribution (dot-dashed), the longitudinal one (dashed) and 
total one (solid line). 



In the Fig. (1) one shows Fz for representative virtualities Q 2 from the latest 
HI Collaboration measurements [12]. The longitudinal and transverse contri- 
butions are shown separately. An effective light quark mass (u, d, s quarks) 
was taken, with the value m q = 300 MeV, and the target radius is considered 
R 2 A = 5 GeV -2 . It should be stressed that this value leads to larger saturation 
corrections rather than using radius ranging over R 2 A ~ 8 — 15 GeV -2 . The 
soft contribution comes from the freezing of the gluon distribution at large 
transverse separation as discussed at the previous section. 

From the plots we verify a good agreement in the normalization, however 
the slope seems quite steep. This fact is due to the modelling for the soft 
contribution and it suggests that a more suitable nonperturbative input should 
be taken. To clarify the role played by the soft nonperturbative contribution 
to F 2 , in the Fig. (2) we plot separately the perturbative contribution and 
parametrize the soft contribution introducing the nonperturbative structure 
function F| oft = C so ft2;~ ' 08 (1 — a;) 10 [13], which is added to the perturbative 
one. The soft piece normalization is C so f t = 0.22. Accordingly, we have used 
just shadowing corrections for the quark sector, taking into account only the 
transverse photon wavefunction and zero quark mass. The integration on the 
transverse separation is taken over 1/Q 2 < r 2 < 1/Qq, with = 0.4 GeV 2 for 
leading order GRV94 gluon distribution. This leads to a residual contribution 
to the soft piece which would come from the transverse separations r 2 < 1/Q 2 . 
It is again verified that the soft contribution is important at small virtualities 
and decreasing as it gets larger. The data description is quite successful. 
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Fig. 2. The Glauber-Mueller prediction for the F% structure function in the rest 



frame. For sake of comparison, one uses quark sector {R\ = 5 GeV 
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and only transverse wavefunction. Radius integration 1/Q < r < 1/Qq an< ^ so ^ 
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Concluding, we have a theoretical estimate, i.e. no fitting procedure, of the 
inclusive structure function F 2 (x,Q 2 ) through the Glauber-Mueller approach 
for the dipole cross section, detecting a non negligible importance of a suitable 
input for the large dipole size region. 



3.2 The structure function Fl 



From QCD theory, the structure function F L has a non-zero value due to the 
gluon radiation, as is encoded in the Altarelli-Martinelli equation (see [14]), 
considering the Breit frame. Experimentally, the determination of the Fl is 
quite limited, providing few data points. Most recently, the HI Collaboration 
has determined the longitudinal structure function through the reduced double 
differential cross section, where the data points were obtained consistently 
with the previous measurements, however being more precise and lying into a 
broader kinematical range [12]. 

In Fig. (3) we present the estimates for the Fl structure function, in represen- 
tative virtualities as a function of £ [5]. For the calculations, it was considered 
light quarks (u, d, s) with effective mass m q = 300 MeV and the target radius 
R\ = 5 GeV -2 . The large r region is considered by the freezing of the gluon 
distribution in this region. Our expression for the observable is then, 
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Fig. 3. The Glauber-Mueller estimates for the Fl structure function. One uses light 
quarks (m q = 300 MeV), target size B? A = 5 GeV -2 and frozen gluon distibution at 
large r. Data from HI Collaboration [12]. 

(I 2 00 1 

F L (X,Q 2 ) =Y^— fd*Tfdz\V L {z,T)\*0™^{x,T>). 
4 7T Z a P m J J 



The behavior is quite consistent with the experimental result, either in shape 
as in normalization. The quantity is less sensitive to the non-perturbative con- 
tent than F2. A better description can be obtained by fine tunning the target 
size or the considered gluon distribution function, however it should be stressed 
that the present prediction is parameter-free and determined using the dipole 
picture taking into account unitarity (saturation) effects in the effective dipole 
cross section [5]. We verify that the rest frame calculation, taking into account 
the dipole degrees of freedom and unitarity effects produces similar conclu- 
sions to those ones using the Breit system. For instance, in a previous work 
[14], the unitarity corrections to the longitudinal structure function were es- 
timated in the laboratory frame considering the Altarelli-Martinelli equation, 
with unitarized expressions for F 2 and xG(x, Q 2 ), obtaining that the expected 
corrections reach to 70 % as ln(l/x) = 15, namely on the kinematical corner 
of the upcoming THERA project. 

The higher twist corrections to the longitudinal structure function have been 
pointed out. For instance, Bartels et al. [15] have calculated numerically the 
twist-four correction obtaining that they are large for Ft and Fl, however hav- 
ing opposite signs. This fact leads to remaining small effects to the inclusive 
structure function by almost complete cancellation between those contribu- 
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tions. The higher twist content is analyzed considering the model [8] as initial 
condition. Concerning F^, it was found that the twist-four correction is large 
and has negative signal, concluding that a leading twist analysis of F L is un- 
reliable for high Q 2 and not too small x. The results are in agreement with 
the simple parametrization for higher twist (HT) studied by the MRST group 

in Ref. [16], where F^ T (x, Q 2 ) = F 2 LT (x,Q 2 )(l + The second term 

would parametrize the higher twist content. In our case, the unitarity correc- 
tions provide an important amount of higher twist content, namely it takes 
into account some of the several graphs determining the twist expansion. 

3.3 The structure function F| c 

In perturbative QCD, the heavy quark production in electron-proton interac- 
tion occurs basicaly through photon-gluon fusion, in which the emitted pho- 
ton interacts with a gluon from the proton generating a quark-antiquark pair. 
Therefore, the heavy quark production allows to determine the gluon distri- 
bution and the amount of unitarity (saturation) effects for the observable. In 
particular, charmed mesons have been measured in deep-inelastic at HERA 
and the corresponding structure function -F 2 cc (x, Q 2 ) is defined from the dif- 
ferential cross section for the cc pair production. 

Experimentally, the measurements of the charm structure function are ob- 
tained by measuring mesons D* ± production [17]. The function F^ c (x : Q 2 ) 
shows an increase with decreasing x at constant values of Q 2 , whereas the rise 
becomes sharper at higher virtualities. The data are consistent with the NLO 
DGLAP calculations. Concerning the ratio R cc = F% c /F2 ) the charm contri- 
bution to F 2 grows steeply as x diminishes. It contributes less than 10% at 
low Q 2 and reaches about 30 % for Q 2 > 120 GeV 2 [17]. 

Once more the color dipole picture will provide a quite simple description for 
the charm structure function in a factorized way. Now, the Glauber-Mueller 
dipole cross section is weighted by the photon wavefunction constituted by a 
cc pair with mass m c . Our expression for the charmed contribution in deep 
inelastic is thus written as 

Ff(x,Q 2 ) =j^ i fd 2 rjdz (|*f?(*,r)| 2 + |*£ M| 2 ) ^f ole (*,r 2 ) 

where l ( z i r )| 2 * s the probability to find in the photon the cc colo dipole 
with the charmed quark carrying fraction z of the photon's light-cone mo- 
mentum with T, L polarizations. For the correspondent wavefunctions, the 
quark mass in Eqs. (2,3) should be substituted by the charm quark mass m c . 
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Fig. 4. The Glauber-Mueller result for the structure function as a function of 
Bjorken variable x at fixed virtualities (in GeV 2 ). One uses charm mass m c = 1.5 
GeV, target size R 2 A = 5 GeV~ 2 and frozen gluon distibution at large r. Data from 
ZEUS Collaboration [17] (statistical errors only). 

Here, we should take care of the connection between the Regge parameter 
x = (W 2 + Q 2 )/(Q 2 + 4 m 2 ,) and the Bjorken variable XBy For calculations 
with the light quarks these variables are equivalent, however for heavier quarks 
the correct relation is XBj = x ( Q 2 /Q 2 + 4 m 2 ), e.g. see [18]:. 

In Fig. (4) we show the estimates for the charm structure function as a function 
of xb] at representative virtualities [5] . In our calculations, it was used charm 
mass m c = 1.5 GeV, target size R\ = 5 GeV -2 and frozen gluon distibution at 
large r. We have verified small soft contribution, decreasing as the virtuality 
rises. There is a slight sensitivity to the value for the charm mass, increas- 
ing the overall normalization as m c diminishes. Such a feature suggests that 
the charm mass is a hard scale suppressing the non-perturbative contribution 
to the corresponding cross section. This conclusion is in agreement with the 
recent BFKL color dipole calculations of Nikolaev-Zoller [18] and those from 
Donnachie-Dosch [19]. 

Regarding the Breit system description, in Ref. [14] we found strong correc- 
tions to the charm structure function, which are larger than those of the F 2 
ones. Considering the ratio R 2 = F2 GM (x,Q 2 )/F2 DGLAF (x,Q 2 ), the correc- 
tions predicted by the Glauber-Mueller approach would reach 62 % at values 
of ln(l/x) ~ 15 (THERA region). Then, an important result is a large devia- 
tion of the standard DGLAP expectations at small x for the ratio R cc = _F 2 CC /F 2 
due to the saturation phenomena (unitarization). With our calculation [5] one 
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verifies that it is obtained a good description of data in both reference systems, 
suggesting a consistent estimation of the unitarity effects for that quantity. 
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